Non-Gaussianity in a Matter Bounce 
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A nonsingular bouncing cosmology in which the scales of interest today exit the Hubble radius in a 
matter-dominated contracting phase yields an alternative to inflation for producing a scale-invariant 
spectrum of adiabatic cosmological fluctuations. In this paper we identify signatures in the non- 
Gaussianities of the fluctuations which are speciflc to this scenario and allow it to be distinguished 
from the results of inflationary models. 
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I. INTRODUCTION 

One of the problems of the inflationary scenario [H, 
the current paradigm of early universe cosmology, is the 
presence of an initial singularity. Such a singularity is 
unavoidable if inflation is obtained from matter scalar 
fields in the context of the Einstein action for space-time 
0. As a consequence, there has been a lot of interest 
in resolving the singularity by means of quantum gravity 
effects or effective field theory techniques. 

A successful resolution of the cosmological singularity 
may lead to a non-singular bouncing cosmology. Such 
non-singular bounces were proposed a long time ago [1]. 
They were studied in models motivated by approaches 
to quantum gravity such as the Pre-Big-Baiig model 14 , 
higher derivative gravity actions (see e.g. 0, 0, 0, [1]), 
brane world scenarios ^9] or loop quantum cosmology 
[lol |. Nonsingular bounces may also emerge from non- 
perturbative superstring theory, as has been investigated 
using tools such as c = 1 matrix theory [TJ] or the 
AdS/CFT correspondence [H, [11]. String gas cosmol- 
ogy Ein approach to string cosmology in which the 
temperature of the universe is non-singular, may also be 
embedded in a bouncing cosmology, as realized e.g. in 
. Finally, non-singular bounces may be studied using 
effective field theory techniques by introducing a curva- 
ture term pi3j | or matter fields violating the key en ergy 
conditions, for exaniple non-conventional fluids EJ, 
or quintom matter jl9|. A specific realization of such a 
quintom bounce occurs in the Lee- Wick cosmology stud- 
ied in H 113] 

In the context of studies of bouncing cosmologies it 
has been realized that fluctuations which are generated 
as quantum vacuum perturbations and exit the Hubble 
radius during a matter-dominated contracting phase lead 
to a scale- invariant spectrum of cosmological fluctuations 



today |22,[2a,[23,|2a| (see also [2g). This yields an alter- 
native to cosmological inflation in explaining the current 
observational data. We will call this scenario the "matter 
bounce". 

It is important to study ways to distinguish the predic- 
tions of the matter bounce from those of the standard in- 
flationary paradigm. One criterion is the tensor to scalar 
ratio which is typically large for a matter bounce [13, [13| . 
In this Letter, we study the distinctive signatures of the 
matter bounce in the non-Gaussianities of the spectrum 
of cosmological fluctuations. 

We find that the amplitude of the non-Gaussianities 
is larger than in simple inflationary models. This is due 
to the fact that the fluctuations grow on super-Hubble 
scales in a contracting universe. This growth also leads to 
a specific signature in the shape of the non-Gaussianities 
which emerges since it is different solutions of the fluctu- 
ation equation which dominate the spectrum compared 
to the case of standard inflation. Both the amplitude and 
shape of the three-point function are predicted indepen- 
dent of any free parameters. 

In the following, we first review the key new features 
which effect the evolution of cosmological fluctuations in 
contracting backgrounds, specifically in the case of the 
matter bounce. Then, we turn to the calculation of the 
non-Gaussianity function f^L and compare our result 
with the predictions of simple slow-roll inflation models. 



II. FLUCTUATIONS IN A CONTRACTING 
UNIVERSE 

To discuss cosmological fluctuations, we write the met- 
ric in the following form (see |27i| for a detailed exposition 
of the theory of cosmological perturbations and [1^ for 
an overview): 



\ri) [(I + 2$)dr/2 - (1 - 2$)dx2 



(1) 
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where 77 is conformal time and $(2;, 77) describes the met- 
ric fluctuations [4^] . 

To follow fluctuations from the contracting to the ex- 
panding phase we will be following most of the literature 
and will use the variable (, which describes the curvature 



2 



fluctuations in co-moving coordinates. If we take matter 
to be a scalar field Lp, then Q is given by 



where 



z 



(2) 



(3) 



a prime indicating a derivative with respect to conformal 
time, and Ti. standing for the Hubble expansion rate in 
conformal time. 

The variable ^ is closely related to the Sasaki- 
Mukhanov variable [1^ [13] in terms of which the ac- 
tion for cosmological perturbations has canonical kinetic 
term: 

V = zC, (4) 
The equation of motion for the Fourier mode Vk of v is 



(k' )vk = 0. 

^ z ' 



(5) 



If the equation of state of the background is time- 
independent, then z ^ a and hence the negative square 
mass term in ([5]) is iJ^. Thus, on length scales smaller 
than the Hubble radius, the solutions of ([5]) are oscillat- 
ing, whereas on larger scales they are frozen in, and their 
amplitude depends on the time evolution of z. 

On super-Hubble scales, the equation of motion ^ for 
in a universe which is contracting or expanding as a 
power p of physical time t, i.e. 



becomes 



which has solutions 



with 



P(2p - 1) _2 



v[j]) ~ 77" 



a = — 
2 



1 1 - 3p 
2I 



V 



(6) 



(7) 



(8) 



(9) 



In the case of an exponentially expanding universe we 
must take the p — s- cx) limit of the above solutions. Hence, 
the solutions for a are a = 3/2 and a = —1/2. The scale 
factor a is proportional to jy^^. Hence 



(10) 



where c\ and C2 are constants. Since 77 — > as i — s- cx), 
it is the constant mode which dominates. This leads 
to the well-known result that the fluctuations of C are 
constant on super-Hubble scales. This also leads to the 



conclusion that the power spectrum of C from quantum 
vacuum perturbations which are initially on sub-Hubble 
scales is scale-invariant since 



Pc(fc,77) 



127r2' 
const , 



(11) 

k^k^^k-^ 



making use of the quantum vacuum normalization and 
the Hubble radius crossing condition. 

In the case of a matter-dominated contraction we have 
ly = —3/2 and hence 



ciV^ + C2T] ^ 



(12) 



where ci and C2 are again constants. The ci mode is 
the mode for which ^ is constant on super-Hubble scales. 
However, in a contracting universe it is the C2 mode which 
dominates and leads to a scale- invariant spectrum [2^[23l. 



Pcik,r^) ^ k^\vkinfa-\rj) 



(13) 



k'\vMm'C-^f 



const . 



^3-1-2 



once again using the scaling of the dominant mode of Vk , 
the Hubble radius crossing condition rjH{k) ~ k~^, and 
the vacuum spectrum at Hubble radius crossing. 



III. NON-GAUSSIANITIES IN THE MATTER 
BOUNCE 

A. Formalism 

In this section we will consider non-Gaussianities in 
the matter bounce model. Specifically, we will focus on 
the amplitude and shape of the three-point function and 
on the function /jvl which is commonly used to describe 
the leading-order non-Gaussianities. 

Non-Gaussianities in single field inflation models were 
first considered in [3l|, [H, [H, [H, IH, and in more de- 
tail in [l^l in the context of single field slow-roll inflation 
models and it was concluded that the non-Gaussianities 
would be small. That larger non-Gaussianities can be 
obtained in multi-field inflation models or DEI inflation 
models was realized in [sl] and [s^, respectively. An 
elegant formalism for calculating non-Gaussianities was 
presented in 40] and extended to the case of generalized 
inflation models in 41]. For a comprehensive review of 
non-Gaussianities the reader is referred to [11] . 

The presence of interactions in the Lagrangian leads 
to non-Gaussianities. We will study them following the 
formalism established in [40] . It is easiest to work in the 
interaction picture, in which the three-point function to 
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leading order in the interaction coupling constant is given 
by 



<c{t,ki)at,k2)c{t,k3) > 



(14) 



= ^ dt' < [C{t,hK{t,k2)C{t,k3),L,r,t{t')] > , 

where ti corresponds to the initial time before which 
there are any non-Gaussianities. The square parenthe- 
ses indicate the commutator, and Lint is the interaction 
Lagrangian (the integral over space of the interaction La- 
grangian density C given below). 

To calculate the non-Gaussianities in the function C, 
we require the Lagrangian density for C to cubic order. 
This has been derived in [401 and takes the form 



£3 = {e'^j)a^CC' + eM{dCr-2e^a^adO(dx) 



,3 xr 



(15) 



where we define x = d and 5 ^ is the inverse Lapla- 
cian. The function / in the last term is 



/(C) 



and e is given by 



H 



(16) 



(17) 



In inflationary cosmology, e is the slow-roll parameter 
and is generally much smaller than order unity, whereas 
in our bouncing cosmology 



3, ,3 
e= -[l + w) ^ - 



(18) 



in the contracting phase [w is the equation of state pa- 
rameter) . 

In the standard cosmological model which describes 
a universe which was expanding starting from an infla- 
tionary phase, C is a constant on scales larger than the 
Hubble radius. On these scales, the spatial derivatives 
are negligible, and thus one sees from ITBl) that the in- 
teraction Lagrangian vanishes. Hence, the integration in 
(fTS]) runs only up to the time of Hubble radius crossing 
in the inflationary phase, i.e. one must just consider the 
non-linear growth of C, inside the Hubble radius. How- 
ever, in a contracting universe C grows on scales larger 
than the Hubble radius. Hence, the integration in (fT5)) 
is dominated by times when the scale is super-Hubble 
(until the bounce time, after which C stops growing on 
super-Hubble scales). As we show below, this leads to 
a difference in the shape of the non-Gaussianities. One 



way to see this is to note that C, is oscillating on scales 
smaller than the Hubble radius, whereas the oscillations 
are frozen out on super-Hubble scales and the growth 
of C, in the contracting phase occurs as the increase in 
the amplitude of a frozen wave perturbation. The very 
different time dependence of C on the scales that dom- 
inate the integral (fTS)) leads to a quite different scaling 
of the integrals, and hence to a different shape of the 
non-Gaussianities. 

There is another important difference between a 
bouncing cosmology and the inflationary model. In sim- 
ple single field slow-roll inflation models the cubic La- 
grangian is suppressed by the slow-roll parameter e ^ 1 
and so the terms proportional to in Eq. psp can 
always be neglected. However, in a bouncing cosmol- 
ogy such as the matter bounce these terms contribute at 
the same order. Moreover, in inflationary cosmology the 
function /(C) is dominated by the first two terms in Eq. 
(dni because on large scales C, is conserved and thus C = 
for the dominant mode. In contrast, in the matter bounce 
it is the last three terms in Eq. (|16p which are dominant 
since the dominant mode of C is growing as 77^"^. Since 
different terms dominate, we will get a different shape of 
the non-Gaussianities. 

The three-point function can be expressed in the fol- 
lowing general form: 



< C(A:i)C(fc2)C(fc3) > 



xA{ki,k2,h) , (19) 



where ki = \ki\ and A is the shape function whose am- 
plitude is (in the case of "local" non-Gaussianities) char- 
acterized by the non-linearity parameter Jnl, where 



(20) 



and Cg{x) is the linear (and hence Gaussian) part of C- 
More generally, in momentum space the amplitude 1^1^^ 
of the non-Gaussianities can be described by 



|6|jV^(fcl,/c2,fc3 



10 ^(fcl,fc2,fc3) 

3 E^kf 



(21) 



As a last preliminary, we note that in the Lagrangian 
formalism, the curvature perturbation variable C. in 
Fourier space can be canonically expressed as 

ari,k) = Uk{v)al + ul{v)a_j: (22) 

with the matter contracting phase mode functions 



Ukiv) = A 



i[l - ik{r] - fjB)] 



exp[ik{r] - fiB)] , (23) 



where fjB is the conformal time when the singularity 
would occur if the matter contracting phase would con- 
tinue to arbitrary densities. The creation and annihi- 
lation operators a and obey the standard canonical 
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commutation relations. The amplitude A is determined 
from the quantum vacuum conditions at Hubble radius 
crossing in the contracting phase. This amplitude de- 
termines the power spectrum of C- If we factor out the 
amplitude and the factor of fc^'^/^, we can define the fol- 
lowing rescaled mode functions: 



(24) 



B. Contributions to Non-Gaussianity in the Matter 
Bounce 

In the following we insert the cubic interaction La- 
grangian ([T5|) into Eq. (fT4|) and calculate the vacuum 
expectation value of the three-point function contributed 
by the interaction terms one by one. We evaluate the 
non-Gaussianity at the bounce time rjB- 

However, before doing that we employ the same trick 
as used in [i^]: In order to cancel the last term in Eq. 
(fT5)) . we make the following field redefinition 



C C-/(C) 



(25) 



Inserting this field redefinition into we find two 

terms: the three point function of the rescaled field ( 
on one hand, and terms in which one factor of C has been 
replaced by /(C) on the other hand. 

To compute the first term, we evaluate the right hand 
side of with an interaction Lagrangian which does 
not contain the last term in p^ . Below, we consider the 
contributions of each remaining term in psp . The second 
term is called the "field redefinition term" . It does not 
involve any integration over time. 

Now we consider the individual terms: 

• Contribution from the Geld redefinition 

Since on large scales the first two terms of /(C) can 
be neglected, we only need to consider the other 
three. Note that this is precisely the opposite of 
what happens in the inflationary paradigm where 
it is the first two terms which dominate. As a conse- 
quence of this difference, the field redefinition term 
leads to a very different shape function. 

Moreover, from the solution of the equation of mo- 
tion for C one can see that there is an approximate 
relation 



c 



(26) 



valid on scales larger than the Hubble radius. 
Therefore, we obtain the following approximate 
form of the redefinition term in momentum space. 



C(fc) ^ C(fc) 



dk 



13 



(27r)3 
9 [k ■%')$■{% 



,3 9 k' -(k- 
|_ ^ 5^ 

^28 fc'2 



k') 



k')) 



k^k'' 



}C{k'K{k-k') ,(27) 



The corresponding shape function is given by (in 
this and the following formulas we keep the factor 
of e explicit since it will allow us to understand 
at which order in e the key differences in shape 
compared to simple slow-roll single field inflation 
models arise) 



— 2^^fcffc| — 2 fcffc|fc^ — kfk^k^y 



e 

32 



(28) 



32nfc| 



fcjfc| + ^ fcf — 2 ^ fcf fc 



^^3 



Now let us turn to the terms which come from in- 
serting the interaction Lagrangian into the right- 
hand side of ((T4)) . These terms all involve an in- 
tegration over time from the initial time until the 
bounce time r\B- They involve a six point func- 
tion of a Gaussian field which yields a cyclic sum 
of products of three two-point functions. Thus, the 
amplitude of the result will be proportional to the 
cube of the two point-function. The Lagrangian to 
be inserted into is the integral over space of 
the Lagrangian density p3|) . Each factor of C in 
(|15p is expanded in plane waves. Making use of 

< C*(fc)C(fc') > = (27r)4fc-3^3(;. _^ k )Pc^(k) , (29) 

we see that the three momentum integrals are ab- 
sorbed by the three delta functions which arise 
when writing the six-point function as a product 
of three two-point functions. The final integra- 
tion over space yields an overall factor (5'^(X]^=i ^i) 
which represents momentum conservation. We will 
demonstrate the steps in computing the six-point 
function for the first contribution, and simply give 
the results in the other cases. 

I Contribution from tlie CC^ term 

The contribution of this term in P^ to the three 
point function is 

{2n)'S{J2h)^^Xlir^B)Xl{VB)Xt,{vB) (30) 



xz 



dvi^ - -)a^Xk,X'k^X'k^ + Sperms + c.c. 



which gives the following contribution to the shape 
function 



A 



^^12 + 24^^^^' 



(31) 
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Note that the shape is different from the contri- 
bution to the shape function of the corresponding 
term in inflationary cosmology, calculated e.g. in 
[4l| . The reason is that in the case of inflation, 
the time integral runs over times during which the 
mode functions are oscillating. Thus, the time in- 
tegral produces a factor of K^^. In our case, the 
integral is over super-Hubble scales and the time 
integration has a very different result. 

• Contribution from the CdCdx term 

A similar calculation shows that the contribution 
of this term to the shape function takes the form 



which we separate first into contributions which arise at 
various orders in e: 



A 



24 



E + i^fp2 { E ^'^^3 - E } 



Once again, the form is different from that of the 
contribution of the same term in inflationary cos- 
mology, for the same reason as explained above. 

Contribution from the Cididjxf' term 

In this case, the contribution to the shape function 
is expressed as 



3 E 



E 

48^ 



E/ ^i^j^k 
i^j^k 



g3 



(33) 



• Secondary Contribution 

One may notice that we have neglected the second 
term of Eq. (fTS]) . Since the form of its shape func- 

tion is approximately taken as kf)(^ h^)' 
contribution from this term is suppressed on large 
scales. 

Finally, summing up all contributions, we obtain the 
following shape function of the three-point correlator, 



A' 



E^; 



3 

i 1 



(34) 



A'" = --yk^ + -yk,k 



32 



(35) 



96 



ffppE^^^i 3^fcf/c 



2j:k!kl 



(36) 



(32) 



96 



5:fcffc| + 3^fcffc: 



which adds up, for our particular value e = 3/2, to 
3 



A, 



256 



j|fc2 1^ E + E ^'i^. 



-66 ^ k^kpl+9 J2 ^tk^^j^. 

i^j=jtk i=itj^k 



(37) 



Comparing with the results for single field slow-roll [40| 
or generalized [4ll | inflationary models, we recognize some 
familar terms (the two last terms in l|37p) and some new 
terms (the terms in the first two lines of ([37)) ). The terms 
which are different from what is obtained in the case of 
inflation arise at second and third order in e. 

Considering the full result ([57)) . we see that the second 
to last term has the largest coefficient and hence dom- 
inates. It is the same term which dominates in simple 
single-field inflation models. Thus, we conclude that the 
dominant term in the non-Gaussianities has local shape, 
and an amplitude which is given and independent of any 
model parameters. The sign is fixed. The new terms 
which are not present in inflationary cosmology are, how- 
ever, not suppressed by more than a factor of order unity. 
Hence, with high quality data they could be seen. 



C. The amplitude parameter 

There are three forms of non-Gaussianity which are 
of particular importance in cosmological observations. 
They are the "local form" , the "equilateral form" and the 
"folded form" , respectively. In single field slow-roll infla- 
tion models, all three are proportional to slow-roll pa- 
rameters and thus are very small. In the matter bounce. 
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the amplitude of the non-Gaussianities is not suppressed 
by slow-roll parameters. Hence, it is clear that matter 
bounces will predict sizable values of these parameters. 

The local form of non-Gaussianity requires that one 
of the three momentum modes exits the Hubble radius 
much earlier than the other two, for example, ki <JC 
k2, k^. Specifically, one is interested in the case when 
the three momentum vectors compose an isoceles trian- 
gle with ki k2 = k^. Then one gets 



\B 



local 
NL 



35 



(38) 



which is negative and of order 0(1). If our predicted 
shape were exactly local (which it is not), then the 



above amplitude would equal the famous f]^'j^^ param- 
eter. Since the matter bounce model predicts a shape 
which is loosely local, one can loosely speaking phrase 
our prediction as 



(•local 

Jnl 



35 



The equilateral form requires fci 



\B 



cquil 

\nl 



255 
~64 



(39) 



In this 



(40) 



The folded form of non-Gaussianity with fci = 2^2 = 
2^3 takes the value 



\B 



folded 
NL 
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(41) 



From the above examples, we see that all of these three 
values of non-Gaussianity are negative and of sizable am- 
plitude. To quantify this statement, we evaluate the re- 
sult numerically setting fc2 = ^3 = 1 and letting fpf^ be 
a function of fci. The physical value of fci runs between 
and 2. 




FIG. 1: The /nl parameter in the Matter Bounce. 




FIG. 2: The shape of non-Gaussianity in the Matter Bounce. 
The vertical axis is the negative of the shape function. This 
figure shows the contribution of the terms of order e. 



D. The Shape of Non-Gaussianity in the Matter 
Bounce 

It is interesting to determine the shape of the non- 
Gaussianities, which has the potential to distinguish dif- 
ferent cosmological models once the data will be suffi- 
ciently accurate. A useful description of the shape is 
given by 



kik2k3 



(42) 



To obtain a better idea of the shape of the non- 
Gaussianities, we have evaluated Eq. (|42)l numerically. 
Our results are plotted in Figures (H) , [31) , 2]) and [S) In 
the figures, we use the following convention: fcs = 1, the 
X-axis is ki, the y-axis is /c2, and the z-axis corresponds 
to the shape A/kik2k3. Figures [2]) , [3]) and[4]) depict the 
shape functions of the contributions to order e, and 




FIG. 3; The shape of non-Gaussianity in the Matter Bounce. 
This figure shows the contribution of the terms of order . 



FIG. 4: The shape of non-Gaussianity in the Matter Bounce. 
This figure shows the contribution of the terms of order e^. 




FIG. 5; The shape of non-Gaussianity in the Matter Bounce. 
This figure shows the contribution of all terms. 

e'^, respectively, Figure ([5|) shows the shape of the total 
contribution. 

For comparison, the shape function of the non- 
Gaussianities in single field slow-roll inflation to leading 
order in the slow-roll parameter e is shown in Figure [S] 
We see that the dominant structure of the shape func- 
tion (modulo sign) is the same in the matter bounce. 
However, it is also clear that the sub-leading correction 
terms in the case of the matter bounce shape function are 
clearly visible. They arise in the terms which are of the 
order and e'^. Therefore, these next-to- leading correc- 
tion terms in the case of slow-roll inflation are suppressed 
by orders of magnitude (namely by e) compared to the 
leading term. 

From the above analysis, we have learned that the am- 
plitude of the non-Gaussianitics of metric perturbations 
predicted in the matter bounce scenario is of order 0(1), 
much larger than in single-field slow-roll inflation mod- 
els. The shape is dominated by a term of local form, 
but there are sizable corrections with which the matter 
bounce model could in principle be distinguished from 



FIG. 6: The shape of non-Gaussianity in single field slow- 
roll inflation to leading order in e. The vertical axis is the 
amplitude of the shape function. 

large classes of inflationary scenarios, including many 
generalized inflation models 41]. 

There are two basic reasons leading to these differ- 
ences. One is that in a bounce model the analog of 
the slow-roll parameter is large, the other is that the 
perturbations outside the Hubble radius are not con- 
served which provides a new origin to generate non- 
Gaussianities. 



E. Squeezed Limit 

It is also interesting to consider the behavior of shape 
function in the squeezed limit when ki = k2 = k and 
fcg ^ [49| 

In this limit, the leading terms of the above three shape 
functions for the contributions to orders e, and e"^ are 
all proportional to 

^ - S ■ (43) 

However, when we sum the three contributions, we find 
that the leading terms cancel and the total shape function 
takes the form 

21 . 

I squeezed = -— /c . (44) 



IV. DISCUSSION AND CONCLUSIONS 

We have calculated the amplitude and shape of the 
non-Gaussianities in the matter bounce model as quanti- 
fied by the three-point correlation function of C. Since in 
this model the fiuctuations grow on super-Hubble scales 
during the contracting phase, different terms in the in- 
teraction Lagrangian dominate the contribution to the 
non-Gaussianities. In addition, there is no slow-roll sup- 
pression of the amplitude. Hence, both the amplitude 
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and the shape are different from what is predicted in in- 
flationary models. 

The ampiitude of the non-Gaussianity is of the order 
0(1). Both the amplitude and the sign are fixed, inde- 
pendent of any model parameters, different from the re- 
sult obtained in the standard paradigm. The amplitude 
which we predict is smaller than the current upper limits 
from the WMAP results [i^, but they are in the range 
of what the Planck satellite will be able to observe |44| 
[50| . The shape function contains contributions which 
are different in shape compared to what is obtained in 
inflationary models, although the dominant term is of lo- 
cal form. These distinctive contributions arise at order 

and in the parameter which in the case of single 
field slow-roll inflation is the slow-roll parameter which 
is much smaller than 1, but which in our case is 3/2. 
Hence, in our case the extra terms in the shape functions 
are not suppressed by more than a factor of order unity 
compared to the contributions which are of similar shape 
to those which dominate in the inflationary paradigm. 

Thus, we conclude that the amplitude and the shape 
of non-Gaussianities of the three-point function provide 
distinctive signatures for cosmological models based on a 
matter bounce. Since the differences in the shape of the 
non-Gaussianities are mainly due to the growth of the 
curvature fluctuations on super-Hubble scales in the con- 
tracting phase, the shape function of non-Gaussianities 
in other scenarios in which adiabatic fluctuations are pro- 
duced in the contracting phase will be similar to what we 
have obtained here. 

Note that we have considered only the case of a canoni- 
cal kinetic term for the matter fields, and wc have focused 
on the adiabatic mode. Including non-canonical matter 



fields and entropy modes will give further contributions 
to the non-Gaussianities which will make the amplitude 
larger in general, as the inclusion of such effects increases 
the amplitude of non-Gaussianities in the case of inflation 
models. 

The formalism we have developed here can be used to 
compute the non-Gaussianities in other co llap sing sce- 
narios, for example the Ekpyrotic model [4^. There 
has been a lot of recent work on non-Gaussian fluctua- 
tions in the multi-field variant of the scenario |46|| . How- 
ever, the adiabatic mode will also contribute to the non- 
Gaussianities, and this contribution can be computed 
easily using our methods. 
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